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Pseudogap phenomena of high-T c cuprates are examined. In terms of AFM (antiferromag- 
netic) and dSC (d I .2_ H 2-wave superconducting) auxiliary fields introduced to integrate out the 
fermions, the effective action for 2D electron systems with AFM and dSC fluctuations is con- 
sidered. By the self-consistent renormalization (SCR), the NMR relaxation rate Tj -1 , the spin 
correlation length f> and the pairing correlation length £<j are calculated. From this calculation, 
a mechanism of the pseudogap formation emerges as the region of dominant d-wave short-range 
order (SRO) over AFM-SRO. When damping for the AFM fluctuation strongly depends on the 
dSC correlation length through the formation of precursor singlets around (tt, 0) and (0, tt) points 
in the momentum space, the pseudogap appears in a region of the normal state characterized 
by decreasing 1/XiT and increasing AFM correlation length with decrease in temperature. This 
reproduces a characteristic feature of the pseudogap phenomena in many underdoped cuprates. 
When the damping becomes insensitive to the dSC correlation length, the pseudogap region 
shrinks as in the overdoped cuprates. 

KEYWORDS: antiferromagnetism, d a .2_ y 2-wave superconductivity, self-consistent renormalization, van-Hove sin- 
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§1. Introduction 

In the underdoped region of high-T c cuprates, spin and charge excitations widely show pseudogap phenom- 
ena in the normal metallic states from temperatures Tpq much higher than the superconducting transition 
temperature T c . Influence of the pseudogap formation shows up in various different experimental probes such 
as NMR relaxation time, the Knight-shift, neutron scattering, tunnenling, photoemission, specific heat, opti- 
cal conductivity, and DC resistivity .Ef In particular, the angle resolved photoemission spectra (ARPESp'tf 
have revealed that the pseudogap starts growing first in the single-particle excitations around (7r,0) and 
(0, 7r) points at T = Tpq. With decrease in temperature, the formation of the pseudogap gradually extends 
in the direction of diagonal in the Brillouin zone such as (tt/2, tt/2) and seems to continuously merge into 
the d x 2_ y 2 gap below T c . From this observation, it is clear that the (tt, 0) and (0,tt) regions in the single- 
particle excitation have a particular importance in the mechanism of the pseudogap formation because the 
pseudogap starts from these regions. In fact (tt, 0) and (0, 7r) regions are known from the ARPES data as 
the region where the quasiparticle dispersion is very flat with strong damping in the underdoped cuprates. 
This seems to suggest that the pseudogap instability is caused because of such flat and damped dispersion. 
We call such (tt, 0) and (0, tt) momenta "flat spot" and the region around them "flat shoal region" and pay 
a particular attention in this paper to understand the mechanism of the pseudogap formation. We note 
that the flat dispersion around the flat spots is also reproduced in numerical caiGulationsQU^ and plays a 
particularly important role in the understanding of the metal-insulator transitiono© 

A remarkable point of the pseudogap is that it appears in both of spin and charge excitations and its 
gap structure in the momentum space is the same as the superconducting d x 2_ y 2 symmetry with continuous 
evolution through T c . This implies that the pseudogap phenomena have close connection to the super- 
conducting fluctuations. Anfitj^j^Qinewhat puzzling point to be stressed is that although the pseudogap 
structure appears in 1/T 1 TJ3BBII3'ILII ) in many cases, I/T2G which measures Ke\(Q = (tt,tt),ui — 0) con- 
tinuously increases with the decrease in temperature with no indication of the pseudogap. This implies that 
the antiferromagnetic excitations are suppressed below Tpq but not in a simple fashion. A related obser- 
vation is the so called resonance peak in the neutron scattering experiments.t3 A resonance peak above 
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the pseudogap frequency sharply grows below T c or Tpq. This peak frequency w* decreases with decreasing 
doping concentration, which suggests that the resonance peak may well continuously evolve into the Bragg 
peak in the AFM phase. This strongly supports the idea that the AFM fluctuations are suppressed around 
U! = but well retained or even enhanced at a nonzero frequency below Tpq. These observations lead to 
the necessity to consider the superconducting and antiferromagnetic fluctuations on an equal footing and 
examine their interplay in detail. 

From theoretical point of view, the origin of the pseudogap is still controvertiaLtir It is important to 
understand why the pseudogap region is so wide in the underdoped region of the high-T c cuprates. In the 
weak-coupling superconductors, pairing fluctuations are usually observed only very close to T c and the BCS 
mean-field description is rather adequate. However, in the high-T c cuprates, the coherence length is only 
several lattice spacings even within a CuC>2 plane, and much smaller than the BCS superconductors. Such 
strong coupling nature of the pairing makes the fluctuation region wider. In addition, the quantum degener- 
acy temperature (coherence temperature T co h ) approaches zero if the Mott insulating state is approached in 
a continuous fashion. This necessarily leads to the vanishing T c . In the underdoped region, such suppression 
of coherence becomes stronger and it could become smaller than the energy scale of the pairing force(~ Tpq). 
With decreasing doping concentration, if the pairing force remains finite, the separation of T c and Tpq is a 
natural consequence. _ _ _ _ _ 

Sometimes such suppressed T c far below Tpc is modelled by the Bose condensation.olij'llj'EZrEfP However, 
this strong coupling aspect cannot be well described by a simple Bose condensation, because the supercon- 
ducting phase appears near the strongly correlated insulator, the Mott insulator. Because of its correlated 
nature which we try to clarify in this paper, the pairing fluctuations do not appear as simple bosons but as 
a complicated interplay with underlying metal-insulator transition, antiferromagnetic fluctuations, d-wave 
symmetry of pairing, and strong momentum dependence of single-particle excitations persisted in the pseu- 
dogap region. An important point often ignored in the literature is the vertex correction, which we discuss 
its effects in this paper. 

In this paper, we make effort to understand the origin of the pseudo-gap behavior in AFM spin excitations. 
We consider the problem under the condition that both of the dSC and AFM fluctuations play roles in 
low energy scales. The pseudogap seen in the spin excitations is consistently understood from the dSC 
fluctuations and precursor effects to the d-wave superconductors. 

Onenlaau renormalization group treatments of electrons around the van-Hove singular points at the Fermi 
leveOEJEJ) show the appearance of attractive vertices for the dSC channel as well as for the AFM, though 
it is difficult to discuss properties at their strong coupling fixed point from them. Learning from this insight, 
in § |2[ we derive the low-energy effective action for the relevant AFM and dSC fields from the electrons 
around the flat spots. In § |^, mode-mode coupling terms are considered in the effective action. We apply 
and extend the SQR (self-consistent renormalization) theory developed by Moriya and his coworkers to 
improve the RPAEJ In § ||, we give numerical results and their physical interpretations, and discuss the 
relation to the pseudogap in high-T c cuprates. We find two characteristic temperatures, Tpq and when 
the dSC order appears at T = 0. The spin pseudogap with suppression in l/T\T appears below Tpq. 
Below this temperature, the dSC correlations become dominant over the AFM correlations and the resultant 
decrease in the paramagnon damping takes place. The AFM correlations are still increasing with decrease 
in temperature in this region. Below the antiferromagnetic correlation starts decreasing while the pairing 
correlations grow as in renormalized classical regime and the superconducting transition occurs shortly. If 
the paramagnon damping is determined from the scattering from the flat shoal region, as expected in the 
underdoped region, Tpq is substantially higher than T*. In the region T* < T < Tpq, the precursor to 
the superconducting state develops the pseudogap structure. When the paramagnon damping is dominantly 
controlled from the region other than the flat shoal region, Tpq merges into . These reproduce the basic 
features of the experimental results. § |^ is the summary of this paper and the problems left for the future 
are also dicussed. 

§2. Low-Energy Effective Action 

Under the condition that AFM fluctuations are the only important ones in the long-length scale of the sys- 
tem, several theoretical and approximate descriptions were developed. For example, in the Mott insulating 
systems, Chakravarty, Halperin and Nelson developed a renormalizatioH. .group scheme for the low-energy 
AFM exciations by employing the d+ 1-chmensional nonlinear-er modelSES Another quite different treat- 
ment of the AFM fluctuations in the literature is the Gaussian approximation of the paramagnetic metallic 
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phase by starting from the Ginzburg-Landau- Wilson action in the overdamped regime due to the coupling 
to the Stoner excitations.E3t3'r2j) 

There exist, however, more complicated cases where other type of fluctuations also come into play in 
the low energy scale. In this case, interplay and competition between different fluctuations seriously affect 
physical properties in general with possibilities of new critical phenomena. The high-T c cuprates provide such 
examples, where pairing fluctuations are expected to compete with AFM fluctuations, because AFM and 
dSC long-ranged ordered phases are observed close each other in the two-dimensional T—S (temperature-hole 
concentration) phase diagram. 

Even when one starts from the path integral with the auxiliary fields describing only the AFM fluctuations, 
in principle, all the other fluctuations are of course also correctly taken into account if the functional integral 
is performed honestly up to infinite order. In princple, it can be done by taking account all the fluctuations 
around the AFM saddle point. It is, however, practically not easy to handle. Physically meaningful results 
are rather obtained by taking account not only of the AFM saddle point but also additional contributions 
from other relevant saddle points separately. This procedure is justified if the overlap between fluctuations 
around different saddle point is properly considered to keep away from double counting of the degrees of 
freedom. To discuss interplay between dSC and AFM fluctuations, here, we treat the AFM and dSC saddle 
points and fluctuations around them on equal footing. 

We consider a 2D strongly correlated electron system and assume that AFM and pairing fluctuations are 
both strong as in the high T c cuprates. The AFM fluctuations are easily taken into account in the standard 
functional-integral technique by the Stratonovich-Hubbard transformation for the local Coulomb repulsion, 
in which integrals over auxiliary fields for the AFM order parameter appear. The Stratonovich-Hubbard 
field for the AFM fluctuations is given by SU(2) symmetric variable <p a . At the moment the microscopic 
process which drives the d-wave pairing flucutations is not competely understood. Although we discuss 
below possible microscopic origins of the pairing saddle point, we first introduce U(l) symmetric auxiliary 
field (\>d on phenomenologcal grounds. Then the partition function represented by the functional integral 
over both of the AFM and dSC auxiliary fields, 4> a and cj> d , respectively reads 

Z = JvaVaV(j) a V^ d V<j) d e^ s (2.1) 
S = - fly] d 2 k(iQ, m - e k )a a (itt m , k)a a (ifl m , k) 

m •* 

+(3^2 d2 1 2y / \T a \TS{iuj n ,q)-(t) l7 (-iuj n ,~q) + i\J\T d \T (w n , q)(f>d(iu n , Q) + d{iuj ni q)cj> d (iuj n , <?)} 
n •* 

+(t> a {iuj n ,q)-4> (T (-iuj n ,-q) + (p d {iuj n ,q)(f) d {iuj ni q) , (2.2) 

where a and a are Grassman fields for electrons, cf> a is the three-component vector field corresponding to the 
spin, and (f> d and 4> d are the pairing fields creating and annihilating a pair of electrons, respectively. (3 is the 
inverse temperature and w n and f2 m are bosonic and fermionic Matsubara frequencies. The bare dispersion 
is described by e^.. Here, we take the nearest-neighbor transfer t and the next- nearest-neighbor transfer 
t', which leads to = — 2t(cosk x + cosk y ) — At' (cosk x cosk y + 1) — //, where /i is the chemical potential 
measured from the flat spots. Further we have introduced the following operators corresponding to spin, 
and annihilating and creating a d x 2_ y 2-w&ve electron pair, 

S(iu) n , q) = / d 2 kd^(ift m - iu> n , k - q)cr ^a v (iVL mi k) (2.3) 
d(iu n , q)=y^l d 2 ka ll {iVt m , fe)o-^ i/ a r/ (-iO m + iuj n , ~k + q)g(k) (2.4) 

d(iuj n ,q) = ^ d 2 kd fM (in m - iut n ,k - q)a v ^d^(-iil m , ~k), (2.5) 

where cr* is the i-th component of the Pauli matrix, and g(k) — cosk x — cosk y . The term containing r CT 
is obtained from the Coulomb interaction term by the Stratonovich-Hubbard transformation for the AFM 
auxiliary fields <p a while the term containing T d represents coupling to the pairing auxiliary field 4> d . 
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To understand the possible origin of the term containing the dSC fluctuations (j)d and (j>d, the one-loop 
level analyses provide a useful insight. On this level, in 2D systems, we first have to consider the contribution 
from the quasi-particles around (tt, 0) momentum, because the (tt, 0) point contribution generates the most 
divergent logarithmic terms due to the Umklapp scattering and also due to the van-Hove singularity. We 
consider two-dimensional interacting electron systems with the Fermi level being located near the flat spots, 
namely (tt, 0) and (0, tt) points. As we mentioned in §1, these flat spots are not only the original van-How: 
singular point but also the points around which flat dispersions are generically observed experimentallyo'Ej 
as well as numerically even when the bare band structure does not suggest the van-Hove singularity at 
these spotsfl-EllIif One-loop renormalization group theories show that at low energies, both of the dSC and 
AFM susceptibilities becpmadiyergent. In fact the dSC channel is most strongly divergent when the nesting 
conditon is not satisfiedJlj'EJ'ElP It implies that the irreducible vertices for the dSC and AFM channels are 
negative at energies low enough compared with the ultraviolet cutoff. Although the one-loop calculations 
are justified only in the weak coupling region and it only suggests that the fluctuations are scaled to strong 
coupling, it is at least clear that we have to consider the two fluctuations explicitly in the competing region. 

After integrating out the fermions in (2.2), the following effective action is obtained, 

S = + S ? + & + S% + S$ + S% (2.6) 
S£ ] = PY,! d 2 qXa\^ n ,q)<j> a {iuj n ,q)-<t> a {-iLa n ,-q) (2.7) 

n J 

S d ] = PYl d 2 qXd 1 (iun,q)Miu n ,q)<t>d( i Un,q) (2.8) 

n J 

(iu> ni ,qi)-(f> a (iuJ n2 ,q-i)<$>a (iu n3 , q 3 ) ■ 4> a (iw„ 4 , q 4 ) (2.9) 
s dd = P u ddY j d2( li d2 <l2 / ^93^(^,91)^(^3, g 2 )^ d! (ia; n 3 ) g3)^ d (ia; n4 ,g 4 ) (2.10) 

ni,fl3,7l3 

S fd = 2 P U <rd Y j rf2? l / d2q 2 j d2 ^34>a( i ^n 1 ,q 1 )-<t>a(i^n 2 ,q2)M-i^n 3 ,~q 3 )M i ^n i ,q 4 ), (2-11) 
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where 



Xcr(iu n ,q) = (C~ 2 + (q-Q? +lM/cl + (Lu n /c a )^ \ (2.12) 

Xd{iu n , q) = (^ 0)_2 + q 2 + ld\u n \/(? d + K/c rf ) 2 ) 1 , (2.13) 

are AFM and dSC dynamical susceptibilities, q 4 = — q x — q 2 — q 3 , Ua = —ri\ — n 2 — n 3 , and Q is the AFM 
ordering wave vector. In this paper we restrict ourselves to the AFM ordering vector at the commensurate 
value (tt, tt), and neglect possible long-range features of Coulomb interaction which may lead to gapful dSC 
excitations instead of the Goldstone mode even in the dSC ordered state, as in the s-wave SC state. 
and £ d are the spin correlation length and the coherence length at the mean-field level, respectively, and 
the phase excitations (Higgs bosons) are not treated separately from the amplitude modes. The velocity of 
spin and pairing collective modes are denoted by c a and Cd- The damping constants are given by 7^ and 7^. 

We first discuss how £ CT and £ d in (2.12)-(2.13) are derived. If we employ the RPA and T-matrix approx- 
imation, i.e., the Gaussian approximation for all the above fields, £ CT and £ d depend on the temperature 
T and the chemical potential /1 measured from the flat spots in the following way: 

& ral m ^? log ^F log 7 — ru" 2 - 14 

\Jt 2 - At 12 T max{(J,,T} 

C ~1 j-log j ,, Tl log j—^r- (2-15) 

Here E c is a ultraviolet cutoff in the energy scale of the bandwidth and we have included coefficients of 
the double logarithms into T a and Yd- If electrons near the (tt, 0) and (0, 7r) points play dominant role for 
the both ordering fluctuations, in the mean-field level, finite-temperature transitions either to the dSC at 
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T = T c ° , or the AFM at T = T»? are expected as described by 
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E c exp 
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forTV>|/i|,|f| (216) 



for \t'\ > T N > \fi 
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E c exp 



for r c °> > | M | (2 17) 



|r d |iogJs /|A»| 



for T c <0) < | M | 



We next discuss how 7^ and 7^ are determined. In principle, they may be determined microscopically 
from the scattering of the collective modes by quasiparticles. In this paper, however, we consider them from 
phcnomenological arguments. The origin of 7 is mainly from contimnuum of the Stoner excitations and 
the amplitude strongly depends on the density of states of low-energy excitations. This low-energy part 
of damping is usually absent if some kind of long-ranged order appears, and it is also suppressed if the 
correlation length gets longer. When only one type of fluctuations with the correlation length £ exists, a 
plausible dependence for long £ would be 7 = 7 <0) £~ v . 

For the choice of 7, here, two possible cases are of interest. In one case, AFM and dSC fluctuations are 
coupled with quasiparticles excited far from the flat spots (n, 0) and (0, n), and they dominate dampings of 
the above collective fluctuations. In this case we expect (p = 0. In the other case, the main source of the 
damping is ascribed to the quasiparticles around the flat spots. Then we expect <p > 0. 

In §4, we give the calculated results of the pairing correlation length, the spin correlation length and the 
NMR spin-lattice relaxation rate for <p = and 1. When only the AFM spin fluctuations exist, a corssover 
between the dynamical exponent z — 2 at low temperatures and the z = 1 at high temperatures in the 
disordered side near the critical point has been treated by assuming <p = by Sachdev et al£3 Recently, 
Schmalian has discussed the case, ip = 1 by the renormalization group method and has obtained a crossover 
between the behavior with z = 1 at the disordered side near the critical point and the z = 2 behavior far 
from the critical region.© In this paper, we show that in the presence of another strong instability, the 
system shows properties significantly different from the case with only one type of fluctuations. We discuss 
this problem further in §4. 

Since other parameter values c„ and a do not have such singular dependences, in this paper, we neglect 
their dependence on t 1 , /1 and T. For example, we can safely regard the velocities, c a and Cd as constants, 
since the velocities remain finite at the transition point. 

Next we proceed to the perturbative evaluation of quartic terms. The detailed calculations of the coupling 
constants are presented in Appendix [A|. Roughly speaking, the AFM-<iSC coupling constant turns out to 
be positive when < /1 < t — 2t' , < — /z < t + 2t' , or T > are satisfied. The other couplings are 
always positive. Udd is the most divergent at the lowest energy. It diverges as In E c / Max{T, fi} when T — > 0. 
This leads to nothing but the breakdown of the expansion of S e f / in terms of the bare Green function for 
small enough \fi\. We, then, introduce phenomenological constant values for the mode-mode couplings as 
mentioned in §4. In short, our theory displayed below corresponds to a treatment of a phenomenological 
spin-d-wave pairing model. _ _ 

Note that 7r-triplet fluctuations can be suppressed when the SO(5) symmetryE3E$ is absent. In fact, 
recently Furukawa et al. have argued within the g-ological one-loop renormalization group treatment that 
for the (71-, 0) contributions at the Fermi level, this is the case even when the strengths of the divergences,.!^ 
AFM and dSC fluctuations are the same, i.e., when the next-nearest neighbor hopping t 1 is equal to zero.EH 
Therefore we assume that the suppressed 7r-triplet fluctuations may be ignored. 

§3. Self-Consistent Renormalization 

Next using the effective action obtained above, we start renormalization process for the mode-mode cou- 
pling terms. This is a similar procedure to the SCR theory developed by Moriya and coworkers. EJ In our 
case the mode-mode coupling terms consist of those between AFM and AFM fluctuations, dSC and dSC 
fluctuations, and AFM and dSC fluctuations. Following the SCR scheme, we have considered the Feynmann 
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Fig. 1. Diagrammatic representation of the present theory. Solid lines denote electronic Green function while thick (thin) long- 
and short-dashed lines are renormalized (unrenormalized) AFM and dSC fluctuations. The diagrams (a) and (b) represent 
the Dyson equations for the AFM and dSC fluctuations, respectively. The diagrams (c) and (d) are the self-energies for AFM 
and dSC fluctuations. We have kept terms up to the one-loop order in the perturbation expansion of them. 



diagrams represented in Fig. |3|. In the usual spin flucutation theories, diagrams containing dSC fluctuations 
are absent. Here we consider the both AFM and dSC fluctuations to understand their interplay. There 
exists some previous attempts related to this line of approach. For example, the coupling of spins to pairing 
fluctuations, was discussed in the calculation of the uniform spin susceptibility near the superconducting 
transitiono' However this approach did not consider AFM fluctuations. It also did not take an explicit 
account of electrons around the (tt, 0) point. 

The calculated results in fact depend on the choice of the parameter values in the effective action Eq. (2.6). 
Among other parameter values in (2.6)-(2.13), only those representing the bare distance from the critical 
point, £j ' ~ 2 and ~ 2 behave singularly at T = even for finite values of fi, and they are still most divergent 
with logarithmic accuracy. Then they crucially determine the qualitative aspects through the dependence 
on the choices of If, /i and T. 

Consequently, the role of the self-consistent renormalization procedure is reduced to the determination of 
£0- and £d from 

C 2 = C~ 2 +J o J coth |f KJmx ff (w, fe) + u ad Im X d(u, k)} (3.1) 

id 2 = iT~ 2 +J a J coth |f krdInix CT (w, fe) + u dd Im Xd (u, fe)] ■ (3.2) 

Our main results obtained numerically will be discussed in the following section. Here, we comment on 
the characteristics of the present SCR theory for the mode-mode couplings. First, as in the SCR theory 
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of spin fluctuations in 2D, the system can be ordered only at T = 0. ft agrees with the Mermin- Wagner 
theoremH for the AFM order. For the dSC, the K-T transition at nonzero temperatures is not possible in 
our formalism in contrast with what should be. This implies that the treatment of the gauge fluctuations 
in the SCR theory is inadequate for the description of the K-T transition at Tkt- However, once all the 
higher energy degrees of freedom are integrated out, the renormalized superfluid stiffness determines the 
temperature scale where the pairing correlation length starts growing strongly. Since Tkt is of the order of 
the stiffness according to the K-T theory .o 1 Tkt should occur close to this crossover temperature in our 
theory, below which the spin correlation length starts decreasing. In our analysis, we take this temperature 
scale as the signature of the K-T transition. 

Another important point is that though the Gaussian approximation for the electrons around the (ttJ^ 
point predicts the coexistence of the AFM and the c?SC, as discussed in § |2[ as well as mean-field arguments,!^! 
the inclusion of vertex corrections due to the coupling between AFM and dSC modes makes their coexistence 
difficult. It is simpliy because the coupling between AFM and dSC modes are repulsive for realistic cases. 
The coexistence appears only at t' = /a = Q. The further conditions for their stability are given by | | < 
\T a \uad/uaa and |r CT | > \Td\u a d/udd for the AFM, and \Td\ > \T a \u a d/u aa and |r ff | < \Td\u a d/udd for the 



dSC, as shown in Appendix B. A quite different type of coexistence of dSC and AFM may in principle 
be possible if different portion of the Fermi surface than (n, 0) and (0, it) contributes to the AFM order 
while <iSC order is realized only through contributions around (7r, 0) part. Since we do not seriously treat 
contributions from (7r/2,7r/2) part, this possibility is beyond the scope of this paper. 

In what follows, we concentrate on the cases where the dSC ground state is realized at T — while 
substantial AFM fluctuations also persist. In our analyses, we take Yd = T a . Then the system shows the 
competing fluctuations of the AFM and the dSC in the normal metallic (but not necessarily Fermi-liquid-like) 
phase below the mean-field transition temperature. 

§4. Results 

Before showing results, we discuss how the phenomenological parameters are determined. For comparison 
with experiments of hole-doped high-T c cuprates shi 



s shawiLbelow, we take the following rough estimates; t 
0.25eV from photoemission data of YBa2Cu307_ K ,Elrt3^ E c — t, c a — 0.5i and r y a — t. Hereafter we 
take the lattice spacing as the length scale and t as the energy unit. Above choices of 7° and c a are 
taken similar to those obtained by fitting experimental data oiJJTiT and T^q of YBa2Cu3C , 7_ :c at high 
temperatures with the staggered susceptibility of the form (2. 12) .c2P At the moment, some uncertainty exists 
in determining the parameter values describing dSC fluctuations. For simplicity and definiteness, we take 
c c i = Cfi and cf^ = 7 ff . We assume that t' and ±l are small compared with those obtained by tight-binding 
band fittingsEj with the photoemission results, H and take comparable values to T c , so that contributions 
from the flat shoal regions are dominant. We note that experimentally observed Fermi surface may be 
reproduced by such parameters due to the strong correlation effectsH For £, a d , we adopt the one- loop 
results (2.14) and (2.15) with proportionality constants determined so that the spin and pairing correlation 
lengths at T = 1 are a unity. Further we take the values of fcrmionic vertices T a = Td which predict T N of 
order of the temperature below which the system enters a critical region. Unfortunately, it is rather difficult 
to reliably estimate the renormalized values of the mode-mode couplings in the available framework. In this 
paper, we leave T a — Yd and the mode-mode couplings as adjustable parameters and pose a question whether 
the experimental results can be totally reproduced within reasonable choices of their values. Particularly, the 
following conditions have to be satisfied at least to make more detailed comparisons with the experimental 
data: First, the onset temperature T„ where correlation length £^ starts increasing rapidly should be close 
to T c in real materials. Second, £ CT at T — T* is similar to that obtained in the experiments in each case. 
Third, the crossover temperature Tpq obtained in the case <p = 1 at which \fT{T starts decreasing gives a 
pseudogap temperature obtained experimentally. 

First, for the case <p = 0, we show the case where the parameter values correspond to YBa2Cu307.LP We 
plot the calculated dSC correlation length 1/T\T and spin correlation length £ a in Fig. 2. As previously 
mentioned, on this level of approximations, the nature of the original Fermi surface determines whether the 
ground state is AFM or dSC, through the bare irreducible particle-hole and particle-particle susceptibilities. 
For finite values of t' or /i, the ground state is dSC in the one-loop results. In this case, the spin correlation 
length first increases down to the temperature T* and then decreases with further decrease in temperature. 
The crossover temperature depends on t', mode- mode couplings and other parameters in the effective 
action. It makes a crossover between the regime T < T, dominated by the dSC renormalized classical 
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Fig. 2. The calculated spin correlation lengths (solid lines), 
correlation lengths (long-dashed lines), We have taken |r a | 
t' = —0.03, which correspond to YBa2Cu307. 



1/TiT (dashed lines) normalized by its value at T = 1, dSC 
= \T d \ = 0.15, u aa = 2.6, u ad = 1.2, u dd = 1.8, fi = 0.03 and 



fluctuations and the thermally fluctuating regime T > T* . 

The spin-lattice relaxation rate T\ of a 63 Cu nuclei, after a proper subtraction of the contribution from 
the uniform part, is evaluated approximately as 



1/TiToc lim I dklmXo-iw, k)/u 



(4.1) 



Here we have not considered the k-dependence of the nuclear form factor seriously, because it does not alter 
the basic feature. From (4.1), we see that \jT\T has basically the same temperature dependence as 
By taking T c comparable to T m , our present result ioiua = suggests that overdoped cuprates which show 
Tpg ~ T* ~ T c as in Tl2Ba2Cu06,r*r HgBa2Cu04+,5c3' are also classified as this class, besides optimally 
doped systems such as YBa2Cu307.Lf We conclude that our calculated results with ip = are totally con- 
sistent with the experimental results of optimally as well as overdoped cuprates. Even in the underdoped 
region, I^-aSr^CuO.jC 3 ! seems to belong to this category. Recent ARPES data of La2_ x Sr a; Cu04 clari- 
fied exceptionally strong damping in the region (±7r/2, ±7r/2) e3 This may well be reproduced by ip = 
class. The absence of clear pseudo gap in the magnetic response in this compound is consistent with this 
interpretation. 

Next we consider ip > cases. Contrary to the case ip — 0, one has to be careful in determining the 
scaling form of the damping rates. If only AFM spin fluctuations exist, it can follow the simple scaling form 
7o- oc When we take ip — 1, the z = 1 scaling holds within the tree level. However, when we treat 

complicated fluctuations containing both AFM and dSC characters tightly coupled each other through S^J, 
the correct theory should involve two different scales £ CT and £<j- Here, note that the lower energy scale (or 
longer length scale) always determines the damping rates of the modes: physically, the strongest fluctuations 
affect the fermionic low-energy spectral weights, which dominate the mode dampings. To take into account 
this feedback effect on the mode dampings, we introduce the following phenomenological relation: 



7<T,d 



(4.2) 



In terms of bosonic excitations, the relaxation times of collective modes should be determined by the time 
necessary to propagate the scale of the longest correlation length, because the damping is not effective as 
far as the excitations are propagating inside such an ordered domain. Thus we take tp = 1 here. 

The calculated £ CT , ^ and 1/TjT are shown in Fig. 3 for the parameter values corresponding to 
YBa2Cii3 06.63- As in ip — cases, with the decrease in temperature, the spin correlation length first 
increases, reaches its maximum value and then decreases due to the quantum fluctutaions introduced by the 
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Fig. 3. The results for the spin correlation lengths (solid lines), 1/T±T (dashed lines), and the dSC correlation lengths (long- 
dashed lines) in the tp = 1 case. We have taken |r CT | = |r,j| = 0.7, u arT = 3.14, = 3.05, u a< i = 1.0, (j, = 0.02 and 
t' = —0.02, which correspond to YBa2Cu3C>6.63- As in Fig. 2, the data of \jT\T are normalized by its value at T = 1. 




renormalized-classical dSC fluctuations. This temperature, which we refer to as T*, is mainly determined 
from t' and the mode-mode coupling constants. On the other hand, 1/T±T behaves in a different manner 
from tp = cases. When the temperature is lowered, £ CT and £d both first gradually increase. In this region, 
1/T\T also gradually increases. With further decrease in temperature, £d starts growing quicker than £ CT and 
then 7 CT starts decreasing much quicker than £ CT at Tpc- This drives spin excitations from overdamped into 
underdamped region, and forms a spin-excitation peak at a finite frequency lu = lu* in S(Q,u>). Namely, the 
spectral weight start transferring from lu = to the peak region around lu* , which makes decrease in 1/T\T. 
A similar crossover was previously obtained ip a numerical calculation near the quantum transition point 
between dSC and AFM ordered phases Jl3'C3c2P At lower temperatures characterized by T„, the dSC fluctua- 
tions go into the renormalized classical regime, which signals the decrease in £ CT . We again interpret as the 
rough estimate of T c . These properties .are also similar to experimental data in underdoped cuprates with 
a pseudo spin gap, such as YBa2Cu40sEr and Bi2Sr2CaiCu208-CP In the present framework, the crossover 
Tpc increases with increasing t' , because it accelerates the dominance of dSC correlations over AFM, which 
induces reduction of 7 CT . Figure 4 indeed shows that Tpo increases with t' . Figure 5 shows the calculated 
results for tp — 1 and fj, = t' = with u aa less than u^, so that £ CT shows a stronger divergence than ^ and 
the ground state is very close to the quantum transition point between AFM and dSC. 

Irrespective of the damping rate exponent tp, low-energy spin fluctuations are in general suppressed due 
to the dSC fluctuations around its ordered state. However, the manner of suppression is different between 
tp = and 1, i.e., between the case that low-energy fermions coupled to collective modes exist and the case 
that those are absent. 

The peak structure in S1Q,U)) around lu* reproduces some qualitative feature in the resonance peak 
observed experimentally .E£|[9 In our treatment, to* is selfconsistently determined from the competition 
between cZSC and AFM and the value lu* is characterized by the dSC gap amplitude. The AFM fluctuations 
are pushed out from the region lower than lu* due to the dSC pseudogap formation. 

Experimentally, the resonance peak appears to be sharpened with the increasing <iSC correlation length ^ 
and the intensity of the peak remarkably increases below the superconducting T c . This implies a rapid growth 
of the AFM correlation length On the other hand, lu* does not strongly depend on the temperature. 
The frequency lu* depends only on the doping concentration. It in fact decreases only when the dpoitig 
concentration is decreased. A similar feature is also seen numerically near the AFM-eiSC transitionca'cZl' 
This feature is beyond the scope of this paper, because, in our framework of (2.12), lu* has to be proportional 
to i;" 1 for small 7 CT . To realize rapidly increasing jcorrelation length with a fixed finite lu*, we need to 
modify the assumed form (2.12) as discussed before^ To reproduce the temperature dependence, it will 
be necessary to correctly take into account dominat incoherent part in addition to the coherent response, 
although our susceptibility (2.12) overlooks such a subtlety. 
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Fig. 4. Spin correlation lengths (solid lines) and 1/TiT (dashed lines) calculated in the case of tp = 1 for \t'\ = 0.02, 0.06 and 
0.1 from top to bottom curves, and dSC correlation lengths (long-dashed lines) from bottom to top curves. We have taken 
|r CT | = |r,j| = 0.7, u a(J = 3.14, u dd = 3.05, u ad = 1.0, /i = 0.02. It suggests that as |t'| become larger, T, and TpQ increase. 




Fig. 5. The results for the spin correlation lengths (solid lines), 1/TiT (dashed lines), and the dSC correlation lengths (long- 
dashed lines) in the ip = 1 case. We have taken |r CT | = \T d \ = 0.7, u aa = 3.0, u dd = 3.1, u ad = 1.0, fi = and t' = 0, so that 
the spin correlation length shows a stronger divergence than the dSC correlation length. 



§5. Conclusion 

The pseudogap region observed in the high-T c is reproduced as the region with enhanced dSC correlations. 
The experimental results for the magnetic excitations are consistently explained from precursor effects for 
the superconductivity. 

From the electrons around the flat spots, (tt,0) and (0, 7r), we have considered the low-energy effective 
action for the auxilary fields for the AFM and dSC. The one-particle irreducible vertices pfaiJppth channels 
have been taken negative values, as suggested by one-loop renormalization group studies.EilliSEP Next, we 
have obtained solutions of the self-consistent one-loop renormalization. The inclusion of vertex correction 
in 2D suppresses finite-temperature transitions. This correctly reproduces the absence of antiferromagnetic 
transition at nonzero T in 2D. The superconducting Kosterlitz-Thouless-Berezinski transition is also absent 
in this approximation in contrast to what should be. This is just an artifact due to neglect of the phase 
modes of the dSC order parameter. In this paper, we take T* as a rough estimate of T c , where is defined as 
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the onset temperature for the renormalized classical regime of the d-wave pairing. Below this temperature, 
spin correlation length starts decreasing. Because of the repulsive coupling between dSG and AFM modes, 
the coexistence of the dSG and the AFM orders does not take place except when both degrees of freedom 
are nearly symmetric as explained in Appendix [b|. 

The pseudogap formation clarified from the interplay of AFM and dSC is summarized as follows: When 
cither the Fermi level \x measured from the flat spot or the next nearest-neighbor hopping t' is nonzero, the 
dSG ground state is realized. Only in a special case of fj, = If = 0, dSC, AFM or coexistent ground state can 
be realized depending on the values of the mode-mode couplings as discussed in Appendix [^. When the dSG 
correlation grows faster, the process of the pseudogap formation depends on the damping exponent tp. If 
tp = 0, namely the damping 7 does not depend on the dSC correlation length, the spin correlation length, £ CT 
and l/( 63 TiT) both reaches its maximum value at T = and then decreases as the temperature decreases. 
On the other hand, when the mode dampings are suppressed near the transition, e.g., tp = 1, \/(T\T) shows 
a faster decrease at Tpq(> T*) while £ CT continues to increase until T*. The pseudogap temperature Tpq 
increases as t' increases and it tends to decrease when the nesting effect of the Fermi surface around the flat 
spot becomes stronger. 

In comparison with experiments in the cuprates, the results for <p = reproduce the feature obtained 
in the optimally doped and overdoped high-T c cuprates, such as YBaaCuaOjLl' and Tl2Ba2Cu06 and 
HgBa2Cu04+5,c3 and even some underdoped cases as I^-^Sr^CuO^EJ For this category, we specu- 
late that quasi-particle excitations around (tt/2,tt/2) point mainly contribute to the finite damping of AFM 
fluctuations. 



On the other hand, in many underdoped high-T c cuprates, 1/{T\T) starts decreasing at Tpc at 
while the spin correlation length measured by T^g shows the Curie- Weiss behavior down to T c .l 




This pseudogap behavior is reproduced by our tp = 1 case, if we take as the superconducting transition 
temperature. In addition, the qualitative similarity between our results for tp = 1 and the experimental 
results suggests that the damping of the AFM and dSC collective modes decreases in the pseudogap regime. 
It means that low-energy fermions around the flat spots mainly contribute to the damping. This is consistent 
with the strong damping of quasiparticle around the flat spot observed experimentally in the underdoped 
region. This flat dispersion is a consequence of the universal character of the metal-insulator transition. 

We clearly need further studies for a more complete understanding of the pseudogap in the high-Tc 
cuprates. After the present work, several important steps still remain. In the-pext step from the present work, 
a more systematic renormalization group treatment would be interesting^ We have concentrated on the 
single-particle excitations only around the flat spots, (tt, 0) and (0, it) assuming their crucial importance in the 
mechanism of the pseudogap formation. However, in the one-loop level, the origin of the flat dispersion and 
strong damping in this momentum region is not fully understood. Experimentally the flatness and damping 
strength appear much more pronounced than the expectation from the one-loop analyses. Numerical analyses 
also support that this remarkable momentum dependence around the flat spots is due to strong correlation 
effects. We have to calculate self-energy corrections as well as the vertex corrections in a selfconsistent 
fashion to clarify the profoundness of such correlation effects. This is clearly the step beyond the one- 
loop level. Near the metal-insulator transition point, vanishing quasiparticle weight has to be led from a 
selfconsistent treatment. This will contribute to clarify how the pairing channel appears and how the flat 
spots are destabilized to the paired singlet. We also note that the dominance of the incoherent weight over 
the quasiparticle weight in the single-particle excitations near the metal-insulator transition may require a 
serious modification in the derivation of the AFM and (ZSC susceptibilities. The Curie- Weiss type form for 
the dynamic spin susceptibility needs to be reconsidered ,EiJ because the spin susceptibility is also determined 
mainly from the incoherent part of the single-particle excitations which we have not considered at all in this 
paper. In addition, we have not considered effects of quasiparticles excitated far from the flat spots. In this 
context, a more comprehensive analyses are required to see the momentum dependence. 
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Appendix A: Calculation of the coupling constants 

Here we give the quartic coupling constants at the lowest energy and the wave vector measured from the 
ordering one. If we take g(k) = sgn(cosA; x — coak y ), after straightforward calculations one obtains 
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1 d 
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(A6) 



£ k Q 

Ucrcr and Udd are always positive. u CT d is positive when < — a < t + 2t', < fx < t — 2i', or T > holds 
approximately. Otherwise it is negative. Consequently, it is positive for typical Fermi surface of high-T c 
cuprates. They are finite at least at finite temperatures, thus, at the energy cutoff. 

Appendix B: Phase diagram at T = for nested cases 

In this appendix, we discuss what type of ground states is established in the presence of positive u a d- We, 
however, regard that u a d is small enough so that u aa Udd > u 2 d . 

To find the lowest excitation spectrum, it is convenient to rewrite (3.1) and (3.2) in the following expres- 
sions, 
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rft. I 



2T 



(B-l) 



When i' 7^ or ^ 7^ 0, within weak coupling arguments (2.14) and (2.15), dSC ordered state is realized at 
T = 0. It is because at sufficiently low temperatures, £j ' ~ 2 + (£^ 2 — £ d ' ' ' ~ 2 )u a d/udd becomes positive, while 
£ d ° ~ 2 + (£~ 2 — ~ 2 )u a d/u C r IJ has a large negative value. For the nested case, u = t' = 0, the insertion of 
(2.14) and (2.15) into (B-l) and (B-2) lead to the followings; 

r 2E c 



Cd 
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"dd 

Merd 
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r d )Alog 2 2 u ad /u dd + (u aa - -24) / — /^^l mXCT ( w ,fc)coth— (B-3) 

-t «dd Jo 7T J (^7T) Z 2i 



r CT )/tiog 2 + C 2M <Td/M ( T ( T + (w d d 



2£ c 



U(j(j J tt J (2tt) 2 2T 



at sufficiently low temperatures. Therefore ground state properties are determined by r CT — TdU^d/udd and 
Td — TuUad/uuu. If the former is positive and the latter is negative, the system has the AFM order. If the 
former is negative and the latter is positive, it has the dSC order. If both values are positive, both orders 
may coexist and thus 7r-triplet fluctuations are also required to consider. If both values are negative, it 
belongs to the quantum-disordered region. 
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